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Abstract: We study kink (domain wall) solutions in a model consisting of two complex scalar 
fields coupled to two independent Abelian gauge fields in a Lagrangian that has U(l) x U(l) 
gauge plus Z2 discrete symmetry. We find consistent solutions such that while the U(l) 
symmetries of the fields are preserved while in their respective vacua, they are broken on the 
domain wall. The gauge field solutions show that the domain wall is sandwiched between 
domains with constant magnetic fields. 
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1. Introduction 

Over the last thirty years or so, the study of solitonic solutions to classical field theories 
has yielded many interesting results of wide relevance to particle physics, cosmology and 
condensed matter physics. The more recent fascination with brane-world models of particle 
physics and cosmology has added new motivation for these kinds of investigations. In this 
paper we will study a simple model of two complex scalar or Higgs fields <pi and <p2 coupling 
to two different U(l) gauge fields and A2^, with the added feature of an exact discrete 
Z2 symmetry under the interchange 1 <-> 2. We will derive solutions to the coupled classical 
field equations that exhibit a kink or domain wall form for the scalar fields. The nature of 
the gauge field configurations self-consistently coupled to the Higgs kinks will be our primary 
object of study. A similar model, without the discrete exchange symmetry was studied some- 
time ago by Witten [1] in the context of a superconducting string solution. The model was 
investigated in more detail by MacKenzie [2] to show that while a symmetry is preserved 
in the vacuum, unexpected topological structures can arise in the interior of a domain wall. 
More recently, Lemperiere and Shellard [3] have reported on the behavior and stability of the 
superconducting currents in Witten's model. 

Our own motivation for this rather abstract investigation lies with the symmetry break- 
ing mechanism proposed in Ref. [4] in the context of brane world models and dubbed as the 
"clash of symmetries" . Briefly, Ref. [4] examines a toy model with Higgs fields in three triplet 
representations of a global SU (3) symmetry, where a discrete permutation symmetry between 
the triplets is enforced. Omitting inessential complications, the vacuum states of the theory 
spontaneously break SU(3) down to SU(2), as well as spontaneously breaking the discrete 
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symmetry. Kink solutions are derived that interpolate between vacua invariant under differ- 
ently embedded SU{2) subgroups. 1 For instance, one can have /-spin asymptotically preserved 
on one side of a domain wall, with V-spin on the other. Although the unbroken subgroups on 
both sides are isomorphic, the different embeddings within the parent group cause additional 
symmetry breakdown at all non-asymptotic points. This additional symmetry breaking is the 
"clash" . The idea is that some of the symmetry breaking we see in our universe might be due 
to such a clash, if our world is indeed a brane in a higher dimensional space. 

This idea is still at the developmental stage; no realistic brane-world model building 
using the clash mechanism has yet been attempted, to our knowledge, though Ref. [6] reports 
on some recent progress. In the course of thinking about the clash of symmetries idea, 
however, an even simpler model field theory with U(l) factors and interchange symmetries 
between the different sectors naturally presented itself as a useful theoretical laboratory. 
The model studied in this paper arose in exactly this way, though, of course, it is also 
entitled to an independent existence as a simple-but-not-too-simple vehicle for the study of 
gauge fields coupled to domain wall Higgs configurations. From this perspective, our work is 
relevant to general studies of superconducting topological solitons, as in Refs. [1-3, 7, 8] for 
example. From the clash of symmetries perspective, the present exercise begins the study of 
the breakdown of local continuous symmetries. 

The rest of this paper is structured as follows: In Sec. 2, the model and the field equations 
are presented. The numerical study of kink solutions to these equations is then presented 
in Sec. 3, while Sec. 4 provides a physical explanation for the solutions. Section 5 contains 
some concluding remarks. 

2. The Model 

Using the the notation of [4] we start with the action for two complex scalar fields 4>i,2 coupled 
to different U(l) gauge fields A\^ 2 . To the overall U(l) x U(l) gauge symmetry we add a Z 2 
discrete symmetry which interchanges the scalars, <p\ <-> <fi 2 and the gauge fields, A± <-> A 2 . 
The discrete symmetry makes the two gauge coupling constants equal in magnitude. The 
Lagrangian is 

C = -Ji^iV - \f^F 2iw + (£>5>i)* (Di^i) + (D%<h)* (D 2 ^ 2 ) ~ Vtfufa), (2.1) 
where 

V{<k,<h) = Al {<f>*<k + ^2 - V 2 ) 2 + \ 2 4>*l<Pl<P*2<f>2. (2.2) 

The covariant derivatives in the Lagrangian are given by 

D lfl = - ieA lfl , D 2il = - ieA 2pi . (2.3) 

1 Qualitatively similar solutions, but to a different theory with a different motivation were discovered by 
Pogosian and Vachaspati in Ref. [5] . 
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The Higgs potential admits two vacuum solutions: 

Vacuum 1 : {ftfa) = v 2 (^0 2 ) = 0, (2.4) 
Vacuum 2: (0*^) = ($$ 2 > = v 2 . (2.5) 

These two vacua are degenerate and are the global minima of the potential for the parameter 
regime 

Ai > and A 2 > 0. (2.6) 

We would like to construct domain wall solutions by requiring the scalar Higgs fields to 
asymptote to different respective vacua on either side of the wall. We will be interested in the 
behavior of the corresponding gauge fields for this kind of Higgs configuration. The boundary 
conditions for the scalars are 

I , / \ i JO z — > — oo , . , , .. J v z — > — oo ,„ 

\<k(z)\ = { and \fa{z)\ =\ , (2.7) 

\ V z — > oo I U z — > oo 

where z is the direction perpendicular to the domain wall. 

It is straightforward to compute the equations of motion for the Higgs fields 

dV 

D afl D%<i> a = ~Qp= "2Ai0a [4>* a 4>a + 4>t<i>b - v 2 ) - A 2 ^ a ^4, (2.8) 

where a, b are either 1, 2 or 2, 1 respectively. The equations of motion for the gauge fields are 
similarly given by 

M fT = 2elm [f a (d» - ieA»M ■ (2.9) 

Since we are going to be looking for static domain wall solutions (i.e. static 1 + 1 
solitons), we search for solutions that depend on z but are independent of all the other spatial 
coordinates and time t. In order to simplify our equations we make use of the temporal gauge, 
Aq = 0. With these choices the equations of motion reduce to 

A u = ^, (2-10) 

e 

A'i x>y = 2e 2 A lx , y Rl (2.11) 

R'l = e 2 {A 2 x + A\ y )Ri + 2X 1 R 1 (Rl + R 2 - v 2 ) + X 2 R 1 Rl (2.12) 

where prime denotes a derivative with respect to z and <j) a = R a {z)e taa ^ z \ The corresponding 
equations for the fields with subscript 2 can be obtained simply by exchanging subscripts 1 
and 2. We see in eqn. 2.10| that the z components of both gauge fields are pure gauge and 
because neither A z {z) nor a(z) couple to the physical degrees of the system, they can be 
neglected. 

The coupled differential equations for this system nominally involves six degrees of free- 
dom (one scalar and two gauge degrees of freedom for each field). However, since the x and y 
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components of each gauge field enter quadratically into their respective Higgs field equations 
of motion, it is possible to rotate to a new basis x and y where one only needs keep track of 
one component of each gauge field. Note that the directions perpendicular to z in which each 
of the gauge fields A\ and A 2 point are independent. We therefore have only four degrees of 
freedom to non-trivially solve for. 

The equations we would like to solve are then 

A" = 2e 2 R 2 1 A 1 , (2.13) 
R'{ = e 2 A\R l + 2X 1 R l {R\ + R 2 2 - v 2 ) + A 2 i?i-R| ; (2.14) 

and 1 *-+ 2. We have suppressed the spatial subscripts on the gauge fields, A. 



For a domain wall solution the scalar fields must obey the boundary conditions in eqn. 2.7. 



Thus, by analyzing eqn. |2.13| we see that the gauge fields are required to have the following 
asymptotic behavior: 

Ai(z-»oo) = e ->/2H*l ^0 and A 2 {z -00) = e-^M _> Q. (2.15) 



We observe that this asymptotic behavior is also consistent with eqn. [HI The values 



of Ai(— 00) and ^2(00) are seemingly unconstrained by any of our differential equations. 
However, note that when z <C —1 for A\{z) or when z S> 1 for A2{z) the solutions become 



linear functions of z, the asymptotic solutions to eqn. 2.13. The linear solutions are due to 
the requirement that R\(z) and R2(z) vanish as z — > —00, +00 respectively (this is because 
we require them to be kink solutions). Thus, the only allowed values of A\{— 00) and ^2(00) 
are either a constant (corresponding to constant asymptotic behaviour) or ±00. Consistent 
with this, we will also impose the boundary conditions 

A' x (z = -00) = const. / and A' 2 (z = +00) = const. ^ 0. (2.16) 

The requirement that these slopes be asymptotically nonzero removes the A\ = A 2 = 



solution from our considerations. Eqn. 2.16 allows the constant slopes for A\ and A 2 to be 
arbitrary. If they are chosen to be unequal, it implies that the corresponding magnetic fields 
B\ and B 2 are unequal, leading to a violation of the symmetry inherent in the problem and 
this may also cause dynamical instability of the brane as will be discussed further in section 
4. Hence, it is natural to choose the slopes to be equal. However, our numerical solutions 
(see Fig. ||) show that even in the asymmetrical situation, slopes of A\ and A 2 are very nearly 
equal. 



The coupled differential equations 2.13 and |2.14j together with the conditions of eqns. 



and§l| constitute our boundary value problem (BVP). 



Since we shall resort to numerics to find solutions it is convenient to transform from 
coordinate z to u which is defined on a compact interval, u G [—1,1], via 



u = tanh(tVAiz). (2.17) 
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With this change of coordinates and the field rescalings 

R a -> vR a , A a -> vA a , 

the equations become 

= 2aR{Ai, 



du 2 
,d 2 Ri 



du 
dRi 



(1 _„ Y ^_ 2 „ (1 _^_ 



aA?i?i + 2i?i(i?f - 1) 



(2.18) 
(2.19) 

(2.20) 



and 1 <-> 2. We have defined q = e 2 /Ai and A = A2/A1. We see that solutions only depend 
on two independent coupling constants and not three. In the case of the pure Higgs model 
with a = (see Ref. [4]), if one takes symmetric (R± + R2) and anti-symmetric (i?i — R2) 
linear combinations of the fields, then the differential equations decouple for the special case 
of A = 4 with analytic solutions, 



Ri = \(l + u) 



R2 



5(1 -«) 



(2.21) 



However, this is not the case in our model for 

We shall also be interested in the energy of the solutions we find, thus we need the 
stress-energy for this system 



T, 



(2.22) 



which for our action yields 

Tfiu = -F llia F£, - F 2fia F? u + 2(D lfl <j> 1 )*(D lu cf> 1 ) + 2{D 2ll <hT (D^fa 



\f^F 1{IV + \f^F 2 , u - (D^<p 1 )*(D lfJ ,(f) 1 ) - {I%<hr(D 2li <h) 



+V(<f>! 



(2.23) 



The energy density is then given by the Too component of the stress-energy tensor. This 
simplifies to 



Too - \ 



{A[(z)) 2 + (A' 2 (z)) 2 ] + (R'^z)) 2 + (R' 2 (z)) 2 



+e 2 A l {z) 2 R 1 (z) 2 + e 2 A 2 {z) 2 R2(z) 2 + V(R U R 2 ), 



(2.24) 



for our static solutions and because of our gauge choice, Aq = 0. Thus, in terms of the 
coordinate u and the rescaled fields the energy density is given by 



-too 
Ait; 4 



(1-u 



2\2 



(d u A l{ u)) 2 + { dMu)? + {duRi{u))2 + {duR2{u)) 2 



4 1 4 

+aA 1 (u) 2 R 1 (u) 2 + aA 2 (u) 2 R 2 (u) 2 + (R^u) 2 + R 2 (u) 2 - l) 2 
+\R 1 (u) 2 R 2 (u) 2 , 



(2.25) 



where X±v sets the scale. 
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3. Numerical Solutions 



The numerical method we employ to solve these coupled differential equations is the 'shooting 
method' using the routines from Numerical Recipes in C++ [9]. One can readily convert our 
system of four coupled second order differential equations to a system of eight coupled first 
order differential equations where the functions are: R\, R 2 , A\, A 2 , R[, R' 2 , A\ and A' 2 . 
This is a boundary value problem with the functions Ri, R2, Ax, A 2 are specified on two 
boundaries but with the functions R[, R' 2 , A^, A 2 not specified on either boundary. The way 
the 'shooting method' works is that one guesses values for the derivative functions at the 
left boundary (u = —1), then with all the functions specified on the left boundary one can 
numerically integrate to the right boundary. One then defines a function which measures how 
well the boundary conditions on the right are matched. Using this goodness of fit function 
one can then use a Newton-Raphson procedure to improve the guess on the left boundary for 
the derivatives. One can then iterate this procedure until the boundary conditions on both 
sides are satisfied to the desired accuracy. One potential difficulty is that if the differential 
equations are reasonably complicated (e.g. non-linear) then the initial guess might need to 
be reasonably good in order for the procedure to converge. 

The differential equations, 2.19] and 2.20| , have poles at u = ±1 when one expresses the 



equations as dX/du = (1 — u 2 )~ 2 x . . . Since we cannot evaluate these equations at u = ±1, 
we set the boundaries at u% = — 1 + e and u 2 = 1 — e. However, because now our boundaries 
are not at u = ±1 (z = ±00) we need to know the asymptotic behavior of our functions 
in order to set up the boundary conditions correctly 2 . For the special case of a = and 



A = 4 the analytic solution, eqn. 2.21 , is known from Ref. [4]. While these are not the correct 
solutions for general a and A, they do exhibit the correct asymptotic behavior as u — ► ±1. 
But as long as e is sufficiently small the correct asymptotic behavior is obtained numerically. 



When we solve our boundary value problem numerically we shall use eqn. 2.21 to set the 
boundary conditions for R\ and R 2 . We also need to know the asymptotic behavior of the 
gauge fields near the boundaries. Substituting A = (l — u 2 )@ into the differential equation for 
A (eqn. we can solve for f3, the scaling behavior in the vicinity of the boundary. Thus 

Ax ~ (1 - u)^ 2 ~ e^ 2 as u -» u 2 = 1 - e (3.1) 
A 2 ~ (1 + n)v / °/2 ^ e v^/2 as u ^ ni = _i + e . (3.2) 

The values of Ai(ui) and A 2 (u 2 ) are not constrained by any of the differential equations and 
are therefore left as free parameters. 

As mentioned before when solving a boundary value problem using the 'shooting method', 
convergence may depend on a reasonably accurate guess of the initial conditions on the left 
boundary. This is the case for our set of differential equations since they have an explicit pole 
at u = ±1. This sensitivity gets worse as e approaches zero. The method we employed to 
address this issue involved starting with a relatively large value of e (e = 0.5) and incrementally 



2 For numerical reasons we can not just set Ri(u\) — 0, Ri(u2) 
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reducing it to its desired value using as the initial guess for the values of the derivatives (R^ 
R' 2 , A'i and A' 2 ) on the left boundary for each step the solution of the previous step. 




Figure 1: Plot of R±, A\, R2, A2 and R\ +R2 against tanh _1 (u) for a = 1, A = 4. The free boundary 
conditions are A\(— 1 + e) = Aa(l — e) = 1 for e = 0.005 which corresponds to left and right boundaries 
at tanfi - (it) = ±3. R\ + R2 is nearly constant for this pair of parameters. 



1.4 




Figure 2: Plot of R±, A\, R2, A2 and R1+R2 against tanh (u) for a = 1, A = 1. Here At (—1 + e) = 
A 2 (l- e) = 1 for e = 0.005. 

In figs. [2] and ^ we see numerical solutions to these differential equations for a variety of 
couplings, a, A and boundary conditions Ai(u\) and A2(«2)- We observe that the gauge fields 
A\ and A2 become linear functions of tanh _1 (n) as u — ► ui and n — > respectively. This 
implies that asymptotically these gauge fields become linear functions of 2, which corresponds 
to a constant magnetic field in the direction perpendicular to both z and x (the direction in 
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Figure 3: Plot of Ri, Ai, R 2 , A 2 and Ri+R 2 against tanh _1 (w) for a = 1, A = 4. Here 1+e) = 1 
and A 2 (\ — e) = 0.5 for e = 0.005. The principal effect of the asymmetrical BC's is to shift the center 
of the brane to the right. 

which the gauge field points), 

By ~ d z Ax(z) = constant. (3-3) 

Thus the asymptotic solution (actually tanh(u) -1 need only be of the order of ±2 to be in the 
asymptotic regime for a typical configuration) on either side of the domain wall is a constant 
magnetic field corresponding to the U(l) fields, which point in uncorrelated directions parallel 
to the domain wall. These solutions have non-zero energy density away from the domain wall 
and thus are infinite energy configurations. The solutions where the magnetic fields are both 
zero corresponds to the choice of a = (i.e. no U(l) gauge fields). 

In figs. [j] and g, we have set A^ui = -0.995) = A 2 (u 2 = 0.995) = 1. With this set of 
symmetric boundary conditions the domain wall is centered at u = 0. In fig. ||]we see that the 
effect of asymmetric BC's is to shift the location of the domain wall. While not apparent in 
the figure the magnitudes of the uniform magnetic field far from either side of the domain wall 
do not exactly match. The choice of e = 0.005 (and e = 0.001 for figs. || and ||) corresponds 
to boundaries at tanh(u) -1 = ±3 (and ±3.8). While e can be made smaller at the expense 
of longer computing time, these values are sufficiently small for our purposes. 

In fig. U we see the energy density of a solution plotted as a function of the transverse 
direction. We see that the energy density is peaked at the center of the domain wall. If we 
treat the asymptotic constant magnetic field on either side of the domain wall as a background, 
then we can compute the energy per unit surface area of the domain wall by subtracting off 
the infinite energy associated with the magnetic field. In fig. [| the surface energy density is 
plotted as a function of A for a variety of values of a (boundary conditions are Ai(—l + e) = 
A 2 (l — e) = 1 at e = 0.001). Observe that this 'renormalized' surface energy density is 
only weakly dependent on the value of the gauge coupling constant a. In fig. || we show the 
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Figure 4: Plot of the energy density against tanh (u) for a = 1 , A = 1 . We have used the boundary 
conditions Ai(-1 + e) = A 2 (l - e) = 1 where e = 0.005. 




01 2345678 

Figure 5: Plot of the 'renormalized' surface energy density against A for a — 0.25, 0.5, 1.0, 2.0 (from 
top to bottom). We have used the boundary conditions Ai(— 1 +e) = A 2 (l — e) = 1 where e = 0.001. 



subtracted energy density corresponding to the constant magnetic field as a function of A and 
a. In both figs. || and || we omit values of A < 1 as they require a significantly smaller value 
for e. 

Our solutions are all plotted in units of tanh(u) -1 and not z since u is the natural variable 
in our system of equations, [2.19| and 2.2C . The length scale tanh(u) -1 is dimensionless and can 
be converted into a physical length by dividing by v^/Xi. The thickness of the domain wall is 
typically ~ (see fig. |j) which can be made arbitrarily small by choosing v\f\\ ^ 1. 



0.025 



0.02 



0.015 



0.01 



0.005 




Figure 6: Plot of the energy density of the uniform magnetic field against A for a = 0.25, 0.5, 1.0, 
2.0 (from bottom to top). We have used the boundary conditions A\{— 1 + e) = A%{1 — e) = 1 where 
e = 0.001. 

4. Discussion 



The numerical solutions displayed above have a natural interpretation in terms of supercon- 
ductivity. Consider, for instance, the currents associated with the U(l) gauge groups, 



J, 



ie[cj>*(d^i) ~ + 2e 2 A iil fa i fa 



where i = 1,2. In terms of the amplitude and phase of fa, the currents are given by 



Ji 



-2eR^d lx a i + 2e 2 A ifl Rl 



(4.1) 



(4.2) 



For our configurations, which depend only on z, and for which eqn. 2.10 holds, it is clear that 
only the x- and ^-components are non- vanishing. They evaluate to 



2e 2 A 



i x,y 



{z)R 2 {z). 



(4.3) 



These steady, z-dependent current densities are uniform supercurrent densities localised to 
the domain wall, with the charged boson fields as the current carriers. 



Equation 4.3 shows that the currents are nonzero only when the gauge field configurations 
are nonzero and vice-versa, so these currents are responsible for dynamically generating the 
magnetic fields. On the side of the wall where Ri ^ 0, the corresponding magnetic field is seen 
to decay exponentially, which is simply a Meissner effect. On the other side of the wall, where 
Ri is tending exponentially quickly to zero, we find the magnetic field Bi tending towards 
a finite, uniform configuration pointing in the plane of the wall. This is consistent with 
the domain wall carrying a uniform sheet of current density pointing in the (0, Ai x , Ai y , 0) 



direction, as per eqn. |4.3| . Our configurations have infinite energy because the domain wall is 
of infinite extent, with current densities uniformly distributed on it. 



- 10 - 



The stability or otherwise of our solutions is an important concern. While a complete 
stability analysis is beyond the scope of this paper, the above considerations suggest that the 
geometrically symmetric solutions such as in figs, [l] and ^ could be stable, whereas asymmetric 
configurations such as those of fig. |3| are not. Let current J\ point in the x-direction in the 
plane of the wall. Then eqn. implies that A\ also points in the same direction, so B\ is 
directed along the y-axis. The Lorentz force on the type 1 charge carriers lies in the negative z 
direction. For sector 2, similar reasoning shows that the corresponding Lorentz force on type 
2 charge carriers points in the positive z direction. For symmetric boundary conditions, these 
forces are equal in magnitude as well as opposite in direction. This is a necessary condition 
for stability. For asymmetric boundary conditions, they are unequal, strongly suggesting that 
such configurations are unstable. 



5. Conclusions 

In order to further explore the idea of the "clash of symmetries" from [4] , we have considered 
a model in which two scalar fields are coupled to their respective gauge fields in a Lagrangian 
which has U(l) x U{\) symmetry. We find consistent static solutions for field configurations 



with the vacuum conditions for the scalar fields specified by eqn. 2.7 and the implied boundary 



conditions for the gauge fields, eqn. |2.15[ We obtain the expected kink-like solutions for the 
scalar fields while the two gauge fields diverge linearly on either side of the domain wall. 

When we consider the idealized configuration of an infinitely thin domain wall, we have 
solutions such that while the U(l) symmetries of the fields are preserved in their respective 
vacua, they are both broken on the domain wall. The gauge fields show that the domain 
wall is sandwiched between domains with constant magnetic fields parallel to the wall. In 
the case of a domain wall of finite thickness, there will be magnetic fields parallel to the wall 
on either side. These are associated with superconducting currents, as in the case of the 
superconducting string solution [1]. 

This model demonstrates that in addition to the breakdown of symmetries on the brane, 
the presence of gauge fields introduces new phenomena, such as the appearance of magnetic 
fields. Background magnetic fields of this kind are reminiscent of the configurations in string 
theory that give rise to non-commutativity of space-time coordinates. It would be very 
interesting to see the logical extension of this model to domain wall solutions with non- 
Abelian gauge fields and to study their dynamical effects in addition to symmetry breaking. 
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